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Ultracompact minihalos have recently been proposed as a new class of dark matter structure. 
These minihalos would be produced by phase transitions in the early Universe or features in the 
inflaton potential, and constitute non-baryonic massive compact halo objects (MACHOs) today. We 
examine the prospect of detecting ultracompact minihalos in gamma-rays if dark matter consists 
of self-annihilating particles. We compute present-day fluxes from minihalos produced in the e'^e~ 
annihilation epoch, and the QCD and electroweak phase transitions in the early Universe. Even 
at a distance of 100 pc, minihalos produced during the e'^e~ epoch should be eminently detectable 
today, either by the Fermi satellite, current Air Cerenkov telescopes, or even in archival EGRET 
data. Within ~1 pc, minihalos formed in the QCD phase transition would have similar predicted 
fluxes to the dwarf spheroidal galaxies targeted by current indirect dark matter searches, so might 
also be detectable by present or upcoming experiments. 
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The identity of dark matter remains one of the key out- 
standing problems in physics. Weakly- interacting mas- 
sive particles (WIMPs) provide a compelling solution [1] 
because their weak-scale masses and cross-sections make 
for a natural explanation of the observed abundance of 
dark matter. As most proposed WIMPs are their own an- 
tiparticles, high WIMP densities would also lead to high 
rates of self-annihilation. Annihilation products might 
then provide indirect evidence of the nature of dark mat- 
ter. Gamma-rays are particularly attractive in this re- 
spect, as they do not suffer the same problems of deflec- 
tion and attenuation as massive, charged species. 

It was proposed [2] that dark matter could be massive 
compact halo objects (MACHOs) of condensed baryons, 
e.g. brown dwarfs or faint stars. These are ruled out as 
the dominant component of dark matter by the cosmic 
microwave background (CMB; [3]), Big Bang Nucleosyn- 
thesis [4] , and microlensing searches [5] . Primordial black 
holes (PBHs) are an alternative, disfavoured by their en- 
ergetic evaporation, gravitational influence [6], and the 
large primordial density perturbations required for their 
production {5 > 30%). For comparison, the initial den- 
sity perturbations from inflation were 5 ~ 10~^. 

Ricotti & Gould [7] proposed a non-baryonic MACHO 
that avoids these constraints, and presents a promising 
new target for microlensing searches. Formation proceeds 
similarly to PBHs, whereby small-scale density perturba- 
tions in the early Universe collapse to a compact body. A 
small-scale power spectrum that is the same as observed 
on large scales [3] provides insufficient power for this to 
occur. Perturbations could however be enhanced by fea- 
tures in the inflaton potential, or phase transitions in the 
early Universe [8] . If a perturbation is small, matter will 
not be sufficiently compressed to form a black hole, leav- 
ing only a compact cloud of gas and dark matter. This 



mechanism requires density contrasts of just 6 > 10^^ 
to proceed, so is far more viable than PBH formation. If 
such ultracompact minihalos (UCMHs) exist they will be 
ultra-dense, and excellent targets for indirect detection of 
WIMPs [9]. 

Here we investigate gamma-ray signals expected from 
UCMHs containing WIMP dark matter. We consider 
UCMHs produced in three phase transitions in the 
early Universe: electroweak symmetry breaking (Tew ~ 
200 GeV), QCD confinement (Tqcd « 200 MeV), and 
e+e^ annihilation {T^e ~ 0.51 MeV). We first discuss 
the masses, density profiles and primordial abundance 
of UCMHs, then WIMP models and annihilation chan- 
nels. We present predicted fiuxes and discuss prospects 
for detection with satellite missions and Air Cerenkov 
telescopes (ACTs). In an appendix, we also give ex- 
plicit predictions from a supersymmetric framework with 
a neutralino WIMP. 

Following matter-radiation equality, ultracompact 
minihalos accrete matter by radial infall [7] as 

M^{z)^5m{^-^y (1) 

where Mh{z) is the total mass of the UCMH at redshift 
z, and Zoq is the redshift of matter-radiation equality. We 
assume that UCMHs at z = grew only until z — 10, 
because by this time structure formation would have pro- 
gressed sufficiently far to prevent further accretion [10]. 
The mass contained within a perturbation at equality is 

5m = f^(^-±^Mn{z^), (2) 

where Mh(2;x) is the horizon mass at the time of phase 
transition X, and = f^cDM/f^m ~ 0.834 [3] is the dark 
matter fraction. We take + 1 = 2.32 x IQ'^Vl^h^ [11], 
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giving Zcq = 3160 with ^y^h"^ = 0.136 from the current 
best fit to the CMB, large-scale structure and Type la 
supernovae [3] . Eq. 2 arises because only the dark matter 
fraction of the horizon mass eventually collapses to form 
a UCMH. At zx, this fraction is simply the ratio of the 
dark matter density to the radiation density, which then 
evolves linearly with the scalef actor. 

During radiation domination, the horizon mass is [12] 
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As T{z) a g,5(z)- 3 oc 5,5(^)^3 (1 + z) [11], with 
R the scalefactor of the Universe and g^s the number of 
effective entropic degrees of freedom, this becomes 
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The horizon mass, temperature and effective entropic 
degrees of freedom at equality can be estimated as 
AfHCToq) = 6.5 X lO^^n^h^)-^ = 3.5 X IQI^Mq [6], 
Tcq = 5.5 n^h^ = 0.75 eV and 5*5 (^cq) = 3.91 [11]. At 
the phase transitions, 5*s(Tew) = 107, g*s{TQGB) ~ 55 
and g*s{Tcc) 10.8 [11], giving (5m{EW.QCD.co} = {5.6 x 
10-19,1.1 x 1O-^O.33}M0. 

The dark matter density profile in an ultracompact 
minihalo is [7] 
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in the radial infall approximation, with the maximum 
extent of the UCMH at redshift z 
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The dark matter in an ultracompact minihalo could 
be further concentrated if baryons collapse and contract 
the gravitational potential. We calculated the density 
profile after adiabatic contraction using the method of 
Blumenthal et al. [13]. This assumes that rM{r) is 
conserved at all r, where M{r) is the mass within ra- 
dius r, and that orbits of the dissipationless WIMPs 
do not cross. We assumed that a fraction F of the to- 
tal halo mass condenses to a constant density baryonic 
core of radius Tcoig- We considered F = 10~^, 10~^ and 
Tcoic/Ryi = 5 x 10~^, 10"'^. The effect of the contraction 
is small for the larger core radius, so we show results only 
for Tcoic/Rh = 10~^. Because the induced contraction at 
r is given by the increase in the baryonic mass within 
r, the contraction caused by a constant density bary- 
onic core is most pronounced around the core's edge. 
This is in contrast to the contraction of halos around 
adiabatically-formed black holes, where the baryons col- 
lapse to a central point, steepening the dark matter den- 
sity profile at all radii. The dark matter density in the 



very centre of a halo does not rise significantly in the 
contraction unless the new baryonic distribution also has 
a pronounced spike at the very centre. 

UCMHs also erode over time as dark matter annihi- 
lates away; being ultracompact and ancient, this effect is 
highly significant. A simple way to estimate the maxi- 
mum density pmax at time t in a halo born at U is [14] 
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where is the WIMP mass and {av) is the annihilation 
cross section (multiplied by the colhsional velocity and 
taken in the zero- velocity limit). We truncate the density 
profiles at r = rcut, setting the density within this radius 
equal to pmax- For UCMHs seen today, t — 13.7 Gyr [3]. 
For non-contracted UCMHs, t-, = i(zeq) = 59Myr [15], 
because they have existed since the time of equality. For 
contracted profiles, — t(10) — 0.49 Gyr [15], as they 
were concentrated at z — 10. 

To estimate the cosmological abundance of UCMHs, 
one integrates the probability distribution of primor- 
dial density perturbations between the UCMH formation 
threshold {5 - lO^^) and the PBH threshold {5 - 0.3). 
We approximate the distribution as Gaussian [16], giving 
a 'bare' relic density (i.e. ignoring post-equality accretion 
and disruption) of 
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Here (t(zx)^ is the variance of perturbations at the time 
of the phase transition. Assuming a scale-independent 
perturbation spectrum of index n, and normalising to 
the perturbations observed in the CMB, a can be ap- 
proximated as [16] 



a(zx) = 9.5 X 10-5 (MH(^x)/10'^g) 



,(l-n)/4 



(9) 



On CMB scales, n ~ 1 [3]. However, the CMB probes 
only a limited number of modes. A different power law 
could plausibly dominate at the small scales relevant to 
UCMH formation; indeed, many inflationary models give 
a running spectral index [6], and phase transitions could 
produce scale-dependent features in the power spectrum 
[8]. The present limit at the scale of PBH/UCMH for- 
mation is n < 1.25 [16]. As they grow by a further factor 
of 290 (Eq. 1) between equality and z = 10, UCMHs 
formed in the e^e~ annihilation epoch could account for 
e.g. 10-"^ of today's dark matter \i n = 1.09 (assuming 
they all survive structure formation). For the QCD and 
electroweak phase transitions, similar abundances could 
be obtained for n — 1.06-1.07. 

The gamma-ray flux from WIMP annihilation, in a 
solid angle Afi and integrated above energy i^th, is 
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where dNf/dE is the differential photon yield from the 
/th annihilation channel. The final integral runs over 
the line of sight to the halo. For a spherically-symmetric 
halo appearing as a point source at distance rf, this is 

(11) 

We use d — 100 pc as our canonical value, but our re- 
sults can be rescaled to any d. With a UCMH mass 
fraction of 10~^ and 3.2 x 10^ Mq of dark matter within 
100 pc of Earth (assuming an NFW halo [17]), we expect 
2 X 10^8 electroweak, 1 x 10^ QCD or 3 e+e" UCMHs 
within 100 pc. At 100 pc, all UCMHs are point sources 
to current experiments. 

In Fig. 1 we show gamma-ray fluxes from UCMHs con- 
taining WIMPs annihilating into either bb or We 
computed these with parton-shower photon yields from 
Pythia 6.4 [18] in DarkSUSY 5.05 [19]. The bb channel is 
common in supersymmetric models, and the fi^ fJ-'~ chan- 
nel is prominent in models which fit the PAMELA and 
Fermi electron excesses [20, 21]. For the bb channel we 
use the canonical cross-section (av) = 3 x 10^^^ cm'^ 
implied by the relic density. For fJ.^ fJ.^ we apply a boost 
factor of 100, corresponding to the minimum Sommer- 
feld enhancement necessary to explain the electron data 
in many models. If a UCMH were situated sufficiently 
nearby however, its compactness might provide the re- 
quired boost factor without needing any Sommerfeld en- 
hancement. 

Despite the increased density, adiabatic contraction 
does not greatly increase the gamma-ray flux. This is be- 
cause the flux profile is dominated by the central region, 
which is not strongly contracted. The Sommerfeld en- 
hancement we used for the ijl^ n~ channel increases rcut, 
making the flux profile less concentrated at the centre 
and therefore more responsive to increases in density near 
rcorc- If fcorc < ^^cut Or Tcore > ^'cut, this cffect is abscut. 

In Fig. 1, we show representative point source sensi- 
tivities of EGRET [22] and Fermi [23] above 100 MeV. 
Fig. 2 gives the expected fluxes as a function of threshold 
energy, allowing for a direct comparison with the sensi- 
tivities of current and upcoming ACTs [24, 25]. 

UCMHs formed in the e+e~ annihilation epoch should 
be observable by either Fermi, MAGIC or HESS, depend- 
ing upon the WIMP mass. They could have already been 
seen by EGRET in some cases, effectively ruling out the 
hb channel up to multi-TeV masses. Given their radial 
flux profiles, UCMHs from the e+e~ epoch as close as 
d ~ 1 pc would even appear as extended sources to Fermi. 
The non-discovery to date of a point source with the spec- 
tral characteristics of annihilating dark matter suggests 
that the amplitude of perturbations generated by e+e^ 
annihilation in the early Universe was likely 5 < 10^^. 
A dedicated analysis of the EGRET and Fermi cata- 
logues (particularly unidentified sources) is required for 
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FIG. 1: Integrated fluxes above 100 MeV for UCMHs annihi- 
lating into either bb or jJ.^ fj,~ pairs at a distance d = 100 pc. 
Curves are shown for different phase transitions and degrees of 
adiabatic contraction. Adiabatically-contracted UCMHs are 
assumed to have a fraction F of their mass collapsed into a 
constant-density baryonic core of radius 10"'^ -Rh- Also shown 
are approximate 5a, power-law, high-latitude, point-source 
sensitivities for 2 weeks of pointed EGRET [22] and one year 
of afl-sky Fermi-hAT [23] observations. SoHd limits indicate 
instruments' nominal energy ranges; see also note 25. 



this statement to be made more definite. Such a study 
might even reveal some UCMH candidates. Limits from 
ACTs are more difficult to obtain, as UCMHs could have 
simply been missed by observing the wrong parts of the 
sky. On the other hand, if microlensing searches detect 
a UCMH from the e+e^ transition, it can potentially be 
followed up by ACTs. 

UCMHs from the QCD phase transition are not yet 
visible at d = 100 pc, but at d = 1 pc their predicted 
fluxes would be comparable to those of dwarf galaxies 
[e.g. 26]. If their abundance and the distance of the 
nearest example from Earth were favourable, they might 
be seen by Fermi or future instruments like the Cerenkov 
Telescope Array (CTA). UCMHs from the electroweak 
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FIG. 2: Fluxes from uncontracted UCMHs at d = 100 pc, as a 
function of the energy threshold of the observing experiment. 
Shaded areas show the regions accessible after a 1 year survey 
by the Fermi-LAT [23], and 50 hr of observation by existing 
and planned Air Cerenkov Telescopes [24]. See also note 25. 



phase transition will not be detectable soon unless some 
lie within the Solar System; in any case, light UCMHs 
would face formation problems from kinetic coupling and 
free-streaming of dark matter. 

These results have important implications. Because 
of Eq. 7, the microlensing profiles of UCMHs containing 
WIMPs could differ from those of Ref. 7. The additional 
annihilation products generated by UCMHs early in their 
lives could have an impact upon the ionisation history of 
the Universe, and photons from the extra annihilation 
might modify the extragalactic gamma-ray background. 
If models explaining the Fermi and PAMELA electron 
excesses are accurate, UCMHs would also inject more 
electrons into the intergalactic medium and increase in- 
verse Compton scattering of the CMB at all wavelengths. 
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Appendix: Flux predictions in the CMSSM 

Fig. 3 shows neutralino annihilation fluxes predicted 
in the Constrained Minimal Supersymnietric Standard 
Model (CMSSM), for UCMHs formed in the e+e" epoch. 
We performed a global CMSSM fit using SuperBayeS [27], 
including CMB constraints on the relic density, accelera- 
tor searches for sparticles and the Higgs boson, the muon 
g — 2, the Bg — Bg mass difference, and limits on rare 
S-decays. Details can be found in Ref. 27. Fluxes show a 
familiar band of high probability from points in the focus 
point region, due to clustering around the canonical an- 
nihilation cross-section compatible with the relic density. 
A lower-probability region is also seen, corresponding to 
models where stau co-annihilation is significant. The en- 
tirety of the allowed CMSSM parameter space should be 
accessible by current instruments if UCMHs were formed 



in the e+e" epoch. Predictions for UCMHs arising in 
the QCD and electroweak transitions look similar, but 
are shifted to lower fluxes as in Figs. 1 and 2. 
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FIG. 3: Expected fluxes at d = 100 pc in the context of the 
CMSSM, for UCMHs formed in the e"^e~ annihilation epoch, 
integrated above 100 MeV. Contours indicate 1 and 2a confi- 
dence intervals. Fits included a range of experimental data, 
and required that the neutralino is the only component of 
dark matter. Predictions from the QCD and electroweak 
transitions look similar, but are ~8.5 and ~17 orders smaller, 
respectively. 



